3D Jet Tomography of Twisted Strongly Coupled Quark Gluon Plasmas 
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The triangular enhancement of the rapidity distribution of hadrons produced in p+A reactions 
relative to p+p is a leading order in (A 1 ^ 3 / log s) violation of longitudinal boost invariance at 
high energies. In A + A reactions this leads to a trapezoidal enhancement of the local rapidity 
density of produced gluons. The local rapidity gradient is proportional to the local participant 
number asymmetry, and leads to an effective rotation in the reaction plane. We propose that three 
dimensional jet tomography , correlating the long range rapidity and azimuthal dependences of the 
nuclear modification factor, RAA(yi,<f),p^\b > 0), can be used to look for this intrinsic longitudinal 
boost violating structure of A + A collisions to image the produced twisted strongly coupled quark 
gluon plasma (sQGP). In addition to dipole and elliptic azimuthal moments of Raa, a significant 
high p± octupole moment is predicted away from midrapidity. The azimuthal angles of maximal 
opacity and hence minima of Raa are rotated away from the normal to the reaction plane by an 
'Octupole Twist" angle, 63(77), at forward rapidities. 

PACS numbers: 12.38.Mh; 24.85.+p; 25.75.-q 



I. INTRODUCTION 

The experimental discovery^ of a new form of 
strongly interacting Quark Gluon Plasma matter (called 
the sQGPpJ) at the Relativistic Heavy Ion Collider 
(RHIC) has raised many new questions about the physics 
of ultra-dense matter produced in relativistic nuclear col- 
lisions. Striking new phenomena observed thus far in- 
clude jet quenching, fine structure of elliptic flow, nuclear 
gluon saturation, and the baryon-meson anomaly. The 
primary experimental parameters available at the RHIC 
are the beam energy, s/s = 20 — 200 AGeV, and the vari- 
ation of nuclear geometry through different B + A com- 
binations. Variations of the global nuclear geometry via 
impact parameter or centrality cuts adds another set of 
important experimental methods and instruments with 
which the geometric configuration of the sQGP can be 
systematically varied. In non central collisions of nuclei, 
the transverse area of the participant interaction region 
has an elliptic shape with an aspect ratio or eccentricity 
that can be varied in a controlled way via global mul- 
tiplicity cuts. In this paper we explore a more subtle 
aspects of the 3D geometry of sQGP produced in non 
central A+A collisions that can be explored by extend- 
ing jet tomography techniques 0] to full 3D taking the 
third "slice" dimension to be the rapidity variable. 

Strong elliptic transverse flow resulting from transverse 
elliptic asymmetry provides the key barometric probe 



of the sQGP equation of state. A common simplifying 
assumption in most hydrodynamic and transport anal- 
yses of the data thus far is that the initial conditions 
are longitudinally boost invariant 4|. However, the ra- 
pidity dependence of elliptic flow data clearly demon- 
strate important deviations of the global collective dy- 
namics from boost invariance away from the midrapidity 
region^, ||J- Even the most advanced 3+ ID hydrody- 
namic calculations 0, M- Q- which relax the ideal 2+1D 
Bjorken initial conditions, have not yet been able to re- 
produce the observed rapidity dependence of elliptic flow. 
The solution may require better understanding of the ge- 
ometry of the initial sQGP in addition to better theoret- 
ical control over the growing dissipative hadronic corona 
at high rapidities^. 

As we show below, the dimensionless parameter that 
controls violations of longitudinal boost invariance in A+ 
A reactions at finite energies is the local relative rapidity 
slope of the low transverse momentum partons 
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where 2Y = log s is the rapidity gap between the projec- 
tile and target. In Au + Au at RHIC ^/s = 200 AGeV, 
6 w 0.45 cannot be ignored. Even in Pb+Pb at the Large 
Hadron Collider (LHC) with ^fs = 5500 AGeV, 8 w 0.28 
will not be much smaller. 

In this paper, we explore some observable consequences 
of such violations of longitudinal boost invariance. We 
predict that there is an intrinsic rotation of the initially 
produced matter in the reaction plane in non-central 
A + A that can induce an octupole moment of the nuclear 
modification factor Raa(v 7^ 0, (f>,p±; b > 0) at moderate 
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to high p± away from mid-rapidities. This effect arises 
from the combination of two basic properties of high en- 
ergy hadron- nucleus reactions with jet quenching dynam- 
ics in A + A. 

One important property is that the rapidity density of 
produced hadrons is observed|lO| to be Triangular rel- 
ative to the density produced in p + p reactions with 
a slope given by ~ S. The second basic property fol- 
lows from the Eikonal-Glauber reaction geometry that 
predicts local variations of the projectile and target par- 
ticipant number density in the transverse plane, Xj_. We 
propose that the long range correlations of the rapidity 
and azimuthal dependence of jet quenching can be used 
to extend current 2D Raa{4>) j e * tomography to "im- 
age" the 3D geometry of the twisted sQGP initial state 
at forward rapidities \rj\ ~ 1 — 3. 

In the next section we begin with a review of the p + A 
triangle and its implications for A + A initial conditions. 
In section IIIII we discuss the role of more realistic dif- 
fuse nuclear geometry on the local participant and bi- 
nary collision densities. In Section II VI we calculate jet 
opacities as a function of 7/, <\>. In Section jet attenua- 
tion through the tilted bulk sQGP density in the reaction 
plane is studied, and in Section IVT1 the azimuthal Fourier 
moments of the nuclear modification factor are computed 
as a function of rapidity. Our proposed new observable, 
the Octupole Twist angle 6*3(77), is computed. We con- 
clude in Section lYlII 

II. THE P+A TRIANGLE 

In p+A reactions, the ratio R p a(ti) = 
(dN p A/dr])/(dNpp/d7]) of produced low p± hadrons 
as a function of rapidity, 77, was found at all energies to 
be "triangular" with the height near the rapidity (-Y) 
of the target nucleus proportional to va oc A 1 / 3 , the 
average number of "wounded target nucleons" [ToL ITl| . 
(Note that we are using the variable 77 to denote the 
particle rapidity rather than the particle pseudorapidity, 
which 77 is conventionally linked to. This is to avoid 
later confusion with the variable y, used to denote the 
transverse vertical coordinate in this paper.) Excess 
nuclear enhancement cx A 1 / 3+<? near the target rapidity 
is due to intra-nuclear cascading of low relative energy 
hadrons produced in the target fragmentation region. 
We will not consider this extra nuclear modification 
further in this paper. 

The projectile proton in the center of rapidity frame 
has rapidity, Y , and interacts at a transverse impact pa- 
rameter, b, on the average with 

v A {h) » a m T A (b) = a m J dz PA (z, b) cc A 1 / 3 (2) 

target nucleons with rapidity —Y. Here (Ji n {\fs) is the 
NN inelastic cross section, and the Glauber profile func- 
tion, T^b) measures the number target nucleons per 
unit area of nucleus with proper density pa{t)- 
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FIG. 1: The approximate triangular form of the ratio of the 
observed D + Av^fy charged particle pseudorapidity distribu- 
tions to the dashed HIJING p + p curve from Fig.l are shown 
for two centrality cuts. The solid curves are predictions us- 
ing HIJING 1.383 Q for D + Au/p + p. The ratio of UA5 
p + p[l3| data to the baseline HIJING p+p predictions are 
shown by filled squares. (In color online) 



The simplest way to understand the origin of the pA 
triangle is in terms of phenomenological string mod- 
els of high energy inelastic interactions. Each wounded 
nucleon[lJJ in both projectile and target is diffractively 
excited into high mass ~ \/s/2 states that decay as 
beam jets into relatively low transverse momentum (j>t ~ 
Aqcd) over a wide rapidity interval. Wounded target nu- 
cleons decay over an interval — Y < 77 < while projectile 
nucleons decay over > 77 < Y. The resulting asym- 
metry of the produced hadron rapidity density dN p A/dr] 
thus grows approximately linearly with va- 

At RHIC 12] , the rapidity triangle was clearly observed 
in the Deuterium+Au control experiment at 200 AGeV 
relative to pp and pp|l3j| as shown in Fig. ^ The ab- 
solute agreement between HIJING calculations and the 
measured multiplicities can be seen in Fig. El the dis- 
agreement at extreme negative rapidities marks the pres- 
ence of intranuclear cascading. Also note that the figures 
show pseudorapidity distributions rather than rapidity 
distributions. These are the same up to a jacobian sup- 
pression factor present in all pseudorapidity distributions 
( 15%) that can be seen in Fig. [21 This jacobian factor 
cancels in the ratio. 

One of the early parton model interpretations of the 
rapidity triangle and its generalization to rapidity "trape- 
zoids" in B+A collisions was given in BGK[HJ. The color 
excitation of wounded baryons proceeds in this picture 
with a parton-parton scattering at a rapidity, y* , dis- 
tributed approximately uniformly between —Y and Y. 
The color excited projectile nucleon then neutralizes by 
producing a string of hadrons from y* < r\ < Y while 
the target nucleon neutralizes by producing a string of 
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FIG. 2: Asymmetric pseudo-rapidity distributions of charged 
hadrons produced in D + Au minimum bias and central 0-10% 
reactions at 200 AGeV from PHOBOS|ll are compared to 
p + p data from UA5|13|. The curves show predictions using 
the HIJING vl.383 code[l^|. (In color online) 




FIG. 3: Schematic illustration of how local trapezoidal nu- 
clear enhancements of the rapidity distributions in the reac- 
tion plane (x,rj,y = 0) twist the bulk initial density about 
the normal in non central A + A collisions, (see eqs. 13181 ) (In 
color online) 



hadrons from —Y < rj < y* . The approximately trian- 
gular rapidity density in this model arises naturally from 
the sum over the va target and one projectile fragments. 

BGK^|| generalized this model to predict "trape- 
zoidal" distributions that would arise when a row of 
vb > 1 projectile nucleons interacts inelastically with a 
row of va > 1 target nucleons. This very basic feature of 
multi-nucleon interactions in pQCD is at the core of the 
phenomenological success of LUND string models such 
as FRITIOF[l3 and HIJING Q. The curves shown in 
Figs.l and 2 are predictions of the HIJING modelyja] us- 
ing version 1.383 that includes the Hulthen distribution 
of nucleons in Deuterium. The magnitude and centrality 
dependence of the charged particle rapidity distributions 
are well predicted by the basic trapezoidal nuclear en- 
hancement of multi-particle production dynamics. We 
note that Wang[l8j studied this rapidity asymmetry also 
as a function of px out to moderately high p± to test 
Cronin versus shadowing dynamics away from midrapid- 
ity. 

While the global features of the p± integrated dN/dij 
are well accounted for, we focus in this paper on the novel 
local implications of the trapezoidal violations of longitu- 
dinal boost invariance for the initial conditions in B + A. 
In a non-central collision at an impact parameter b, the 
local enhancement of bulk particle production relative to 
p+p at a particular transverse coordinate xj_ is approx- 
imately given in the BGK model by 

dRBA M _ dN BA /d V d 2 x ± 

d 2 x± w x ^; b J - dNPP/dll 

«^(x x -b/2)^ + Mxx + h/2 )X±n. . (3 ) 

Even if p + p were longitudinally boost invariant with a 
constant dN pp /drj, the variation of va and vb with xj_ 
leads to local violations of longitudinal boost invariance 
in non-central B + A. The local deviations are controlled 
by 




FIG. 4: As in Fig. 3 from another perspective. The straight 
segments in each rapidity slice are for emphasis only. Real- 
istic diffuse spherical geometries curve these lines (see Fig. 7 
below). (In color online) 



*(x x ;b) = ^=^, (4) 

i.e. the local participant asymmetry vt — vp at x±. 

For very large energies, Y 3> 1, Bjorken longitudi- 
nal boost invariance is restored. However for Au+Au 
at RHIC, 7^5, A 1 / 3 ~ 6, and there is an order unity 
variation of S across the transverse profile of the inter- 
action region that leads to the twisted ribbon geometry 
of the local density distribution as illustrated schemati- 
cally in Figs. 3,4. In these figures, the coordinate normal 
to the reaction plane , y, is set to zero, and the surface 
dRAA is illustrated as a function of x — xt in different rj 
slices. In reality only the fixed x slices have trapezoidal 
shape while the rj slices are rounded due to the nonlinear 
dependence of the Vi on transverse coordinates. 

At the global dN/dr] level, the integration over the 
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transverse coordinates leads to 

dN BA /drj 



Rba(ji\ b) = 



dNPP/drj 



±(N A + N B ) + ^L(N B -N A ) , 



(5) 



where N A = J d 2 x±v A (K±; b) and N B are the total num- 
ber of wounded target and projectile nucleons interact- 
ing at impact parameter b. What is remarkable about 
eas. (|3l5[l and readily seen on in Figs 3,4 is that for sym- 
metric A + A reactions, the global ratio may appear to 
be approximately boost invariant, but the local density 
retains an intrinsic rapidity asymmetry that can be char- 
acterized by a rotation twist in reaction plane 



T 



XT} 



d 2 ( dN 
drjdx °^ \drjd 2 x.± 



(6) 



where x and y are the coordinates in and out of reaction 
plane. A useful property of this characterization is that 
it is independent of the normalization of dN and inde- 
pendent of the multiplicative rapidity modulation factor 
dNPP/dn. 

For a sharp sphere geometry with 



v oc v^R 2 - (x ± b/2) 2 - y 2 



(J) 



the xr\ twist in the center x — y = of the oval interac- 
tion region is independent of r\ and given by 



T xv (x.± = 0) 



1 



2b 



Y 4R 2 - b 2 



(8) 



The twist would of course be absent if nuclei were uni- 
formly thick. The magnitude for realistic spherical nuclei 
depends on the actual diffuse Wood Saxon geometry that 
we use for our numerical estimate below. However, (jSJ 
shows that its typical magnitude of T xn is 1/(RY) in 
minimum bias reactions. 

We also note that the trapezoidal distribution of low 
Pt multi-particle production in D + Au follows more gen- 
erally if the collinear factorized QCD mini-jet dynam- 
ics of BGK is replaced by the kr factorized gluon fu- 
sion mechanism in the Color Glass Condensate (CGC) 
model^, |2£j. Quantitative differences between these 
mechanisms arise at moderate pr > 2 — 5GeV which 
are currently under investigation at RHIC[ll l2l). They 
are connected with different physical approximations to 
Cronin enhancement and gluon shadowing p2^. 

The empirical trapezoidal form © in both BGK and 
CGC cases differs from local equilibrium initial condi- 
tions based on firetube or firestreak models 0, 0, H3- 
While such models also predict an effective tilt in the xrj 
plane, those corresponds to possible initial rapidity shear 
of the hydrodynamic fluid at fixed rj. In our approach, 
there is no initial rapidity shear and the collective longi- 
tudinal "fluid" velocity vanishes in the comoving (y = rf) 
frame over the whole transverse region. However, there 
is an initial gradient of the density along the beam axis 



that can induce a fluid shear during later evolution. Po- 
larization has been proposed as an observable sensitive to 
rapidity shear j24[. We concentrate here on generic trape- 
zoidal (BGK, HIJING, or CGC) initial "twisted sQGP" 
conditions because they follow directly from QCD parton 
dynamics and well constrained by p + A phenomenology. 



III. DIFFUSE NUCLEAR GEOMETRY 

To take the diffuse nuclear geometry into account, we 
recall that the local projectile B and target A participant 
densities is given in Glauber theory by 



dN 



part 



dxdy 
dxdy 



= T B (r+)(l-e 



_ p -<Xi„T A (r_) 



T A (r-)(l-e-' r ^ TB ^) 



(9) 



where r± = y (x ± -|) 2 + y 2 and a w 42 mb at y^i = 

200 AGeV. Rare jets are produced on the other hand 
distributed in the transverse plane according to the the 
binary collision number distribution given by 



dNv 



dxdy 



a H T B (r + )T A (r^) 



(10) 



where (Th(p±,ti) is the pQCD jet cross section. For 
B = A collisions, the binary distribution is obviously 
symmetric in the transverse plane at any impact param- 
eter, unlike the participant distributions. Nuclear shad- 
owing and Cronin effects can break this symmetry, but 
for the present study we assume here that at high enough 
Pt , these initial state nuclear effects can be neglected. 

Jet quenching at a given rapidity in azimuthal direc- 
tion (j) depends on a transverse coordinate line integral 
across the sQGP density profile at that rapidity 3]. The 
initial sQGP density is assumed here to be proportional 
to the local participant density. In the BGK model this 
is then approximately given by 



dN n 



C 



dxdyd-q 2Y 



exp[ 



a 



6(Y-\ V \) 



' " rt (Y-r))+ (Y ■ 



dxdy 



dxdy 



V) 



(11) 



where we have introduced a Gaussian rapidity envelope 
to model dNPP/dr] at RHIC. With Y = 5, a n = 3 
and C w 24, the xj^ integrated dN g /dr] x 1/3 fits the 
BRAHMS central Au + Au — > 7r+ , ir~ ) data[25| as shown 
in Fig. 5. 

Figures El and show the 3D diffuse participant ge- 
ometry of the twisted sQGP initial conditions 1)11(1 in 
the transverse plane in several impact parameter and 
rapidity slices. For b = there is complete azimuthal 
symmetry at all rapidities of course. For finite b the 
initial sQGP density shifts to the left at forward ra- 
pidities and to the right in backward rapidities with 
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FIG. 5: A Gaussian fit of the 7T rapidity density in centrai 
Au + Au at 200 AGeV is shown compared to BRAHMS pflj 
data. (In color online) 
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FIG. 6: Contours of the twisted sQGP initial density in the 
transverse plane (x,y) in different rapidity r\ = —3,0,3 and 
impact parameter b = 0, 6, 9 fm slices. Note the opposite 
transverse shifts at r\ and —r\. (In color online) 




FIG. 7: Transverse profile of the twisted sQGP at b = 6 
Fm as a function of transverse coordinate x at y = for 
r\ — ±2, ±4. The normalization is to the local transverse 
participant density fm~ 2 . The overall decrease of the density 
away from midrapidity is due to the approximate Gaussian 
envelope in Fig. 5. (In color online) 
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our conventions. The shifts are small < 2 fm, but dis- 
tinct. There is also a small increase in the eccentricity, 
e(r?) = ((a; - (x)) 2 ~y 2 )/((x- (x)) 2 +y 2 ), of the twisted 
sQGP density away from midrapidity as shown in Fig|8] 



IV. AZIMUTHAL DEPENDENCE OF THE 
TWISTED SQGP OPACITY 

Unlike the participant scaling bulk sQGP density, the 
binary jet distribution, (|10f> . does not manifest a rapidity 
asymmetry at least at enough high px where shadowing 
is expected to be small. The jet production points are 



FIG. 8: The rapidity dependence of the elliptic eccentricity 
of the twisted sQGP. (In color online) 



therefore not rotated away from the beam direction in 
this limit. This is illustrated in Fig|5] 

Jet quenching due to energy loss in the medium de- 
pends on the opacity and hence density of the plasma in 
the direction of propagation 0, H(|. The mismatch be- 
tween the distribution of jet production points and the 
rotated distribution of the bulk sQGP matter therefore 
induced a peculiar azimuthal dependence of the weighed 
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FIG. 9: A schematic illustration of the orientation of the initial twisted sQGP 111! relative to the symmetric binary collision 
distribution HOH of jet production points projected onto the (a;, rf) reaction plane. The extra azimuthal asymmetry and octupole 
twist of jet quenching results from the relative rotation of these two distributions. (In color online) 



opacity 



Xa(x ,y ,r), 4>i b) = c a / dtt a 

J T 



dxdydrj 



{x +tcas{4>),Vo+tsm{4>),b) . (12) 



Different a = —1,0,1 weights correspond to different 
mechanisms of energy loss. Elastic scattering energy loss 
through longitudinally expanding 4] (or static) sQGP 
matter corresponds to a = — 1 (0). Inelastic radiative 
energy loss through longitudinally Bjorken expanding (or 
static) sQGP matter corresponds, on the other hand, to 
a = (1). The results turned out numerically to be 
within 10% in all cases. Therefore, we drop from now on 
the a index and show results only for a = 0. 

We shown in Figs. EH the 4> dependence of the 
opacity for different b and 77. In all cases the opacity is 
normalized to the maximal xq = yo = b = 77 = opacity 



Xo = Xo(0, 0,0,0) 



(13) 



In Fig. ^3 xo(0, 0, 77, 4>)/xo IS shown for a jet produced 
at the origin x = y = for b = 3, 6, 9 fm. For each b, we 
consider the azimuthal variation at 77 = 0, ±2, ±4. First, 
note that the magnitude of the opacity is proportional 
to the decreasing central density as either b or \tj\ in- 
creases. Second, the amplitude of the elliptic second mo- 
ment (cos(2</>)) increases as the eccentricity of the sQGP 
increases with b as seen in Fig|SJ Third, the maximum 
opacity at rj = is in directions perpendicular to the im- 
pact parameter vector where there is more matter and 
minimum parallel to it. These features and trends are 
the familiar ones in the conventional standard Bjorken 
scenario p^. 

However, the new twists caused by the intrinsic longi- 
tudinal boost violating 5 — 0(1) initial conditions illus- 



1 


b = 3fm 


1 


b = 6fm 


1 






0.9 




0.9 




0.9 








it/2 Jt 3rc/2 




0.8 




0.8 




0.8 




jt/2 it 3jc/2 








0.7 




0.7 


AA 


0.7 


0.6 




0.6 




0.6 


0.5 




0.5 




0.5 


0.4 




0.4 




0.4 


0.3 




0.3 




0.3 


0.2 




0.2 




0.2 


0.1 




0.1 


. "* -. ~ ' 


0.1 
















b = 9fm 

T) = 

ti = -2 

u = 2 

- - ■ T| = -4 

- - ■ T| = 4 




7t/2 71 3rt/2 
$ (Azimuth) 



jt/2 Jt 371/2 
1 (Azimuth) 



Jt/2 jt 3ji/2 
if (Azimuth) 



FIG. 10: Opacity \ as a function of azimuthal angle <f> for 
different values of rapidity and impact parameter. The jet 
is assumed to be born at the origin of the transverse plane. 
The functions have been normalized to xo, the opacity felt 
by a jet propagating along the x-axis after being born at the 
origin of the transverse plane, at midrapidity and with impact 
parameter zero. (In color online) 



trated in Fig. 9 are (1) the emergence of a long range ra- 
pidity asymmetry between +77 and —77 with a 180° phase 
shifted azimuthal dependence and (2) an azimuthal pat- 
tern that has odd harmonics in addition to quadrapole 
t>2 = cos(2<fi) moment. It is clear that there is a sig- 
nificant high pt Vi — cos(cf>) dipole moment away from 
midrapidity. This is manifest in the difference between 
x(0 = 0, 77) and xi^i v) as expected from Fig. 9. 
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FIG. 11: Averaged Opacity \ as a function of azimuthal an- 
gle 4> f° r different values of rapidity and impact parameter. 
The initial position of of the jet seeing the opacity is aver- 
aged over the binary distribution. The functions have been 
normalized to xo, the opacity felt by a jet propagating along 
the x-axis after being born at the origin of the transverse 
plane, at midrapidity and with impact parameter zero. (In 
color online) 



The most interesting new and subtle feature in Fig. 
10 , however, is that the azimuthal angles of maximum 
opacity for rj ^ are no longer at tfi — ir/2 and 37r/2, as 
at rj = 0, but rotated slightly away from the normal to 
the reaction plane. 

The splitting of the maxima is the unique new signa- 
ture of the Bjorken violating initial conditions considered 
here. Both v\ and t>2 at high rapidities are well known 
and expected in nuclear collisions due to directed and el- 
liptic hydrodynamic flow of the bulk matter (see Stoecker 
in ref.2 and |23|). However, the long range rapidity anti- 
correlation of the azimuthal rotation of opacity maxima, 
that as we show below correspond to jet quenching min- 
ima, is the smoking gun that we propose as the test of 
twisted sQGP initial conditions. 

To test the robustness of this signal, we next average 
over all the initial jet production points using the binary 
collision number distribution 



< X >= 



1 



dN B in/dr] 



dx Q dy ( x o, yo, b)x (14) 



The results shown in Figure ^] show all the the same 
qualitative effects as in the center jet results in FigllOl 
but slightly diluted due to the extra averaging over Xo. 

In the case that we consider here, the relatively nor- 
malized dA^Bin actually has no 77 dependence. At mod- 
erate px the interplay between Cronin and shadowing 
could give rise to a rj dependence that we will report in 
a subsequent paper. 



V. AZIMUTHAL DEPENDENCE OF Raa 

The observation of jet quenching in different rj slices 
makes it possible to extend 2D jet tomography to 3D. 
The observable is the nuclear modification factor 



Raa = 



dN AA /drjd 2 p T 
T AA da pp /dr/d 2 p T 



(15) 



This is the ratio between the spectrum of high pt hadrons 
produced in A + A and the inclusive differential cross 
section in proton-proton collisions scaled by the binary 
transverse density, T AA (b) = J d 2 x.T A (r + )T A (r^) = 
N bin (b)/a H . 

In general, the computation of R AA (r), <p,p±)) requires 
averaging over the power law pQCD parton spectra, us- 
ing the distribution of energy loss as in the Gyulassy- 
Levai-Vitev (GLV) or Baier-Dokshitzer-Mueller-Schiff 
(BDMS) theory 0, and integrating over hadron fragmen- 
tation functions. This is beyond the scope of the present 
paper. Since our aim here is mainly to introduce and 
illustrate the qualitative features of twisted sQGP initial 
conditions, the simple phenomenological jet absorption 
model used in ref.[28j will suffice for our purposes. In 
this model the survival of a jet in direction </3 relative to 
the reaction plane is simply assumed to be exp(— kx(4>)- 
Averaging over the production points and taking into ac- 
count the particular geometric dependence of the line 
integral for opacity (|12|) for our geometry, the nuclear 
modification factor is given by 



R AA (ri,4>;b) 



A^Bi: 



d 2 Xo ^Bin (x0)6)e -« X (X ,,^,6) 

d l x 



We simplify the computation of k by simply fitting the 
observedy,] approximately pt independent central colli- 
sion R AA fn 0.2 in Au + Au at y/s = 200 AGeV. (This 
is achieved if k ~ 0.25 in units where the central opacity 
for b = is xo = H// 7 ™ 2 )- 

The computed jet survival probability, Raa, as a func- 
tion of azimuthal angle for different rapidities and impact 
parameters can be seen in Figure^l The solid black line 
in the figure is the R AA = 0.2 independent of 4> is base- 
line fit for 77, b = . The main features to note are the 
same as in Figs. 10,11 but with opacity maxima replaced 
by slightly shifted Raa minima. First, the magnitude of 
Raa increases with both increasing impact parameter as 
well as 77. Next, note that at midrapidity, r\ — 0, there 
is a symmetry between (ft = and <fi = ir as usual. The 
elliptic distribution at b > but r\ = leads to the usual 
minimum of the survival probability at (f> = tt/2, 37t/2. 

However, the reflection symmetry is broken away from 
mid rapidity because the jet source and bulk sQGP 
are shifted relative to each other. As expected from 
Figs. 10, 11, the minimum survival probability is rotated 
away from the reaction plane normal in opposite direc- 
tions in the forward and backward rapidity region. To 
better quantify this effect we introduce in the next sec- 
tion the octupole twist observable. 
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FIG. 12: The jet survival probability, Raa, as a function of 
azimuthal angle (j> for different values of rapidity and impact 
parameter. Note the shifted minima of Raa for 77 = ±2 near 
= f ' T" ano - ^ ne overa U f phases shift characteristic of the 
twisted sQGP initial conditions as in Figs. 10,11. (In color 
online) 



VI. THE OCTUPOLE TWIST AND FOURIER 
DECOMPOSITION 



One way to visualize the azimuthal asymmetry of Raa 
is to compare the jet survival probability in forward and 
backward rapidity on the same polar plot as shown in 
Figure 1131 For each fixed rapidity slice and impact pa- 
rameter, we normalize R AA {ili & b) by its minimum Raa" 
over <p and plot the resulting Raa/RaT as a function of 

For b = collisions, this relatively normalized ratio is 
a unit circle independent of <j> and 77. It is shown as a 
thin circle. For b > 0, but 77 — the transverse density 
profile has an elliptic shape centered at x = y — with its 
longest axis oriented normal to the reaction plane as seen 
in Fig. 9. This leads to the highest survival probability in 
the reaction plane (f> = 0, 7r, as shown by the thick solid 
black line in Figll3l 

For both b > and 77 > 0, the bulk matter is shifted 
toward ip = n (in our conventions) leading to a smaller 
survival probability than in the <fi = direction as shown 
by the solid (blue shifted) pear shape curve in the figure. 
The long range twisted sQGP initial conditions lead on 
the other hand for T) = — \rj\ < to an azimuthal pattern 
(red shifted dashed) that is simply the reflected version 
of the 77 > curve about the y axis. The pair of minima 
of both curves are indicated by the arrows. 

We define an "octupole twist" angle, #3(77,6) to help 
quantify the magnitude of these long range rapidity cor- 




FIG. 13: A polar plot of the relatively normalized Raa/RaT 
for 6 = 6 comparing two 77 = ±2 slices (blue solid, dashed red) 
to the unit circle corresponding to b = (thin black). For b = 
6 but 77 = (thick black), the minima are at </j> = 5, t? and 
maxima are at <j> — 0, n. However, for both b > and \r]\ > 0, 
twisted sQGP initial conditions lead to an octupole "pear" 
deformation with the minima shifted to d> = 



83 



3tt _ 

2 2 ' 

The octupole twist angle, #3(77), is a measure the long range 
in rapidity correlations caused by intrinsic local 0(A 1 ^ 3 /Y) 
violations of longitudinal boost invariance in A + A. (In color 
online) 



relations as 

3 M) = 4>{R7t{ri,b)) - c/>(RT A n (-V,b)) (17) 

The minima in the forward region are tilted to (it — 63) /2 
and (3-7T + 63) /2 while in the backward region they are 
tilted in the opposite direction. 

Note that the reaction plane orientation can be deter- 
mined unambiguously experimentally via standard low 
Pt bulk directed and elliptic flow v\ , «2 @ . Jet tomogra- 
phy of high pt hadrons provides a complementary short 
wavelength probe of the reaction dynamics and can also 
be decomposed into azimuthal harmonics, as 



RA A (r), < t ) iP±;b) = a(l + 2}^v n cos(n<p)) 



(18) 



The numerical evaluation of the first four moments, 
v n (rj,b), are shown in Figures 1151 and 1161 In addition 
to the usual v%, v-i and even V4 moments, the new feature 
shown here is the existence of an octupole V3 — (cos(30)) 
harmonic that is responsible for the pear shape deforma- 
tion of Raa in Fig- 13. The magnitude is relatively small 
though, on the order of 7J4. However, V4 for low px has 
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FIG. 14: The octupole twist angle 6-i(rj) of the minimum jet 
survival direction as in Fig. 13 for different values of the 
impact parameter are shown. (In color online) 



FIG. 16: a(r),b) versus r\ for different values of the impact 
parameter. (In color online) 



b = 3fm 



b = 6fm 



b = 9fm 




-4 -2 2 4 ' -4 -2 2 4 ' -4 -2 2 4 
Rapidity (r,) Rapidity (r,) Rapidity (r,) 



FIG. 15: Plots of the Fourier moments of Raa against rapid- 
ity for different values of the impact parameter. Note that 
the third and fourth moments are shown multiplied by 3. (In 
color online) 



already been successfully measured 0, and with a large 
luminosity upgrade and advanced detectors, a measure- 
ment of U3 at n = ±2 may thus become possible at a 
future RHIC IlJ^. 

The predicted peak of the high pt \vi \ at r\ ~ 2, (with 
opposite sign relative to the spectator v±) is another char- 
acteristic feature that can be used to test twisted sQGP 
initial conditions via 3D jet tomography. 



VII. CONCLUSIONS 

In this paper we applied the local trapezoidal multi- 
particle production BGK model to predict the magni- 
tude of intrinsic longitudinal boost violating features of 
the bulk matter produced in A + A collisions. We noted 
that this is a generic feature of both local mini-jet based 
dynamics as encoded in LUND string and HIJING mod- 
els as well as factorized production models such as 
the CGC. The resulting geometrical picture is that the 
sQGP initial conditions has the expected transverse ellip- 
tical shape for non central collisions, but the bulk density 
is twisted away from the beam axis in the reaction (x, r/) 
plane. 

In contrast, the binary collision distribution of hard 
pQCD jet sources remains reflection symmetric in the 
transverse plane oriented parallel to the beam axis at 
least in the approximation that Cronin and shadowing 
can be neglected at high enough px- These approxima- 
tions will be relaxed in a subsequent paper to explore 
further how the px dependence 3D jet tomography of 
twisted sQGP initial conditions can be used to further 
test different models of gluon saturation. 

The generic displacement of the jet sources and bulk 
matter away from central rapidities was predicted to pro- 
duce a novel jet absorption pattern that deviates from 
what is now known at mid-rapidities. The most intrigu- 
ing is the predicted "octupole twist" of the minimum of 
Raa away from the normal of the reaction plane and 
a peak of the jet V\ near rj ~ ±2. These can be further 
quantified in terms of a non- vanishing octupole harmonic 
of the azimuthal distribution. Even though the predicted 
i>3 is small, on the order of V4 it may be measurable at 

RHIC II m. 

The long range in rapidity Arj ~ 4 — 6 correlations in- 



10 



duced by twisted sQGP initial conditions are expected to 
persist even at LHC energies. This is because the dimen- 
sionless parameter, S ~ A 1 / 3 / log(s), that controls this 
phenomenon depends on the logarithm of the collision 
energy. Logarithms are very slowly varying functions of 
the energy and thus the rate of reduction in 8 with energy 
will be concomitantly slow. 
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